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Charge Order with Structural Distortion in Organic Conductors : Comparison 
between #-(ET) 2 RbZn(SCN) 4 and a-(ET) 2 I 3 
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Charge ordering with structural distortion in quasi- two- dimensional organic conductors 0- 
(ET) 2 RbZn(SCN) 4 (ET=BEDT-TTF) and a-(ET) 2 I 3 is investigated theoretically. By using 
the Hartree-Fock approximation for an extended Hubbard model which includes both on-site 
and intersite Coulomb interactions together with Peierls-type electron- lattice couplings, we ex- 
amine the role of lattice degrees of freedom on charge order. It is found that the experimentally 
observed, horizontal charge order is stabilized by lattice distortion in both compounds. In par- 
ticular, the lattice effect is crucial to the realization of the charge order in #-(ET)2RbZn(SCN)4, 
while the peculiar band structure whose symmetry is lower than that of #-(ET)2RbZn(SCN)4 in 
the metallic phase is also an important factor in a-(ET)2l3 together with the lattice distortion. 
For a-(ET)2l3, we obtain a phase transition from a charge-disproportionated metallic phase to 
the horizontal charge order with lattice modulations, which is consistent with the latest X-ray 
experimental result. 

KEYWORDS: charge order, organic conductor, extended Hubbard model, electron-lattice coupling, 
Hartree-Fock approximation 



1. Introduction 

The family of quasi-two dimensional organic conduc- 
tors (ET) 2 X (ET=BEDT-TTF) exhibit a variety of in- 
teresting physical properties such as superconductivity, 
magnetism and charge order (CO). 1,2 (ET^X consist of 
stacking layers of donor ET molecules and monovalent 
X anions. CO phenomena are observed in several mem- 
bers of (ET^X, which have a 3 / 4-filled 7r-electron band, 
and have been studied intensively since they have an im- 
portant role to our understanding of electron correlation 
effects in low-dimensional systems. 

6>-(ET) 2 RbZn(SCN) 4 3,4 and a-(ET) 2 I 3 5 ' 6 are typi- 
cal compounds that are known to exhibit CO. 6- 
(ET)2RbZn(SCN)4 shows a first-order metal-insulator 
transition at T c = 200K and a spin gap behavior at 
low temperatures much below T C J The transition is ac- 
companied by lattice distortion which changes its struc- 
ture from the #-type in the metallic phase to the so- 
called #d-type in the insulating phase. The CO forma- 
tion below T c has been directly observed by NMR exper- 
iments. 3,4 Several experiments 8 12 such as Raman scat- 
tering 10 and X-ray scattering 11, 12 measurements show 
that the horizontal-type CO is formed in this compound. 
In the metallic phase, X-ray experiments indicate a short- 
range CO with long periodicity, which is different from 
the horizontal-stripe state at low temperatures. 11, 12 The 
NMR measurement also shows slow fluctuations of CO 
in the metallic state. 13 Similar charge fluctuations are 
observed in #-(ET) 2 CsZn(SCN)4, which suggests the co- 
existence of charge modulations with different wave vec- 
tors 14, 15 although no long-range CO has been observed 
in this compound. 

On the other hand, a-(ET) 2 l3 shows a metal-insulator 
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transition at T = 135K, 16 and a spin gap behavior 
is observed below T c by the magnetic susceptibility 
measurement. 17 Charge ordering in a-(ET)2l3 has been 
suggested theoretically by the Hartree-Fock approxima- 
tion 18,19 for the model that takes account of the full 
anisotropy of transfer integrals in this compound. Exper- 
imentally, the CO among ET molecules is confirmed by 
NMR experiments 5,6 and Raman spectroscopy. 20 More- 
over, it has been observed that the charge dispropor- 
tionation exists even in the metallic state. 20,21 The 
horizontal- type CO has been recently and directly ob- 
served by the X-ray scattering study, 22 which also shows 
that the transition is accompanied by a structural dis- 
tortion. 23 25 

Although the compounds #-(ET) 2 RbZn(SCN)4 and a- 
(ET) 2 l3 have COs with very similar spatial patterns at 
low temperatures, the natures of the CO transitions are 
quite different. In #-(ET) 2 RbZn(SCN)4, the transition is 
of first order and accompanied with a large lattice dis- 
tortion and large discontinuity. On the other hand, the 
lattice distortion in a-(ET) 2 l3 is relatively small and, 
although the transition is of first order, the hysteresis 
of the specific heat is substantially smaller than that of 
<9-(ET) 2 RbZn(SCN) 4 . 26,27 In fact, the recently observed 
photoinduced melting of CO in these compounds 28 shows 
a clear difference in the photoinduced dynamics, which 
is considered to originate from different roles of electron- 
lattice couplings in COs. 

Theoretically, CO phenomena have been investigated 
by using the extended Hubbard model including on-site 
(U) and intersite (Vij) Coulomb interactions. 29 41 The 
stability of various CO patterns in (ET) 2 X has been dis- 
cussed first within the Hartree approximation by consid- 
ering the realistic band structures of (ET) 2 X. 2,29 The 
horizontal-stripe CO is shown to be stabilized in 0- 
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(ET) 2 RbZn(SCN) 4 and also in a-(ET) 2 I 3 . In particular, 
the horizontal CO is shown to be more stable in the Od- 
type structure than in the #-type one, which suggests 
the lattice effects have an important role to realize the 
horizontal CO in (9-(ET) 2 RbZn(SCN) 4 . For <9-(ET) 2 X, 
COs with long periodicity are also considered by sev- 
eral authors 32 ' 34 40 and they discussed their relevance to 
anomalous charge fluctuations in the metallic phase. 

In the previous paper, 40 we have studied the effects of 
lattice distortion on CO in #-(ET) 2 X within the Hartree- 
Fock approximation by taking account of three Peierls- 
type electron-lattice couplings explicitly, which are de- 
duced from the X-ray scattering experiment. 42 The re- 
sults show that the horizontal CO is stabilized by the 
lattice distortion, which results in the structural change 
from the #-type to the ^-type. This is consistent with 
the experiments on #-(ET) 2 RbZn(SCN)4 and also con- 
sistent with the recent numerical exact-diagonalization 
study. 41 

In this paper, we discuss the role of each electron- 
lattice coupling to the formation of CO in detail for 9- 
(ET) 2 RbZn(SCN)4. Furthermore, we consider the lattice 
effects in a-(ET) 2 l3 and compare the results with those 
for (9-(ET) 2 RbZn(SCN) 4 . It is found that although the 
lattice degrees of freedom are crucial to the realization 
of horizontal CO in (9-(ET) 2 RbZn(SCN) 4 , the effects are 
relatively small in a-(ET) 2 l3 where the band structure 
and Coulomb interactions have an important role. At 
finite temperatures, we compare the free energies of var- 
ious CO patterns and investigate the CO phase transi- 
tions in both compounds. This paper is organized as fol- 
lows. In §2, the extended Hubbard model and the Peierls- 
type electron-lattice couplings for #-(ET) 2 RbZn(SCN)4 
and a-(ET) 2 l3 are introduced. We show the results of the 
Hartree-Fock calculations at zero and finite temperatures 
and discuss their relevance to the experimental results in 
§3. §4 is devoted to the summary and conclusion. 

2. Model and Method 

2. 1 Formulation 

We consider the extended Hubbard model written as, 



(1) 



+ uJ2 n n n n + Y v ^o n i n 3 + Y ^Y u h ' 

* (ij) (ij) 

where (ij) represents the summation over the pairs of 
neighboring sites, cj^c^) denotes the creation (annihi- 
lation) operator for an electron with spin a at the ith 
site, rii a = c\ a Ci a , and rii = + n^. The electron 
density is fixed at 3/4 filling. The electron-lattice cou- 
pling constant, lattice displacement and elastic constant 
are denoted by a^j, Uij and K iy j, respectively. For the 
intersite Coulomb interactions V^j, we consider nearest 
neighbor interactions V c for the vertical direction and 
V p for the diagonal direction as shown in Fig. 1(a). We 
further rewrite the parameters for the lattice degrees of 
freedom by introducing new variables as yij = aijUij 
and Si j — afj/Kij. The actual form of y^j together 



with the transfer integrals Uj for #-(ET) 2 RbZn(SCN) 4 
and a-(ET) 2 l3 are given below. We apply the Hartree- 
Fock approximation 



( C ia C j^ , ) C ja /C ' 



'(cLc ia ) (2) 



to eq. (1) and the obtained Hamiltonian is diagonal- 
ized in fc-space by assuming unit cells of various mean- 
field order parameters. We consider four types of or- 
der parameters according to the alignment of hole-rich 
molecules, namely, 3-fold, diagonal, horizontal and verti- 
cal COs which are schematically shown in Fig. 2. As for 
the spin degrees of freedom, we use three spin configura- 
tions in each stripe-type CO which are identical to those 
of ref. 29, while the antiferromagnetic state between hole- 
rich and -poor sites is considered in the 3-fold CO. The 
ground-state energy is calculated by solving the mean- 
field equation self-consistently together with the lattice 
displacements, which are determined by the Hellmann- 
Feynman theorem. The energy per site is given by 



Y V ij( n i)( n j) + V ij( C L C j^)(c] a Ci a ) 
(ij) (ij)° 



(3) 



+£ — 

as) 2Slj 

where Ei\^ a and rip is a band index, an energy eigen- 
value of the mean-field Hamiltonian and the Fermi dis- 
tribution function, respectively. N is the total number of 
sites. We notice that the presence or absence of the Fock 
terms in eq. (2) does not alter the qualitative results in- 
cluding the lattice displacements although the energy of 
each self-consistent solution is lowered by the Fock terms. 

At finite temperatures we calculate the free energy 
within the Hartree-Fock approximation. The free energy 
per site is written as 

F = 1 (jiN tot -^Y ln ( 1 + e*P{-P(E llur - /i)}) 



Zkcr 



(4) 



(i3)v 



(ij) 



where /i, N to t and (3 is the chemical potential, total num- 
ber of electrons and inverse of temperature, respectively. 

2.2 0-(ET) 2 RbZn(SCN) 4 

The structure of <9-(ET) 2 RbZn(SCN) 4 in the metal- 
lic phase and that in the CO phase are shown in Figs. 
1(a) and 1(b), respectively. In the metallic phase, there 
are two kinds of transfer integrals t p and t c for diagonal 
and vertical bonds, while the CO transition doubles the 
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Fig. 2. Spatial patterns of COs considered in the Hartree-Fock 
calculations. (a)3-fold, (b)diagonal-stripe, (c)horizontal-stripe, 
and (d)vertical-stripe COs. The hole-rich and hole-poor sites are 
denoted by the solid and open ellipses, respectively. 



connected bonds. These three kinds of electron-lattice 
couplings cause the modulations yij of transfer integrals 
that are experimentally observed. For simplicity, they are 
assumed to be independent and we do not consider any 
other electron-lattice coupling. Then, the transfer inte- 
grals in the distorted structure are given by, 



Fig. 1. Schematic representations of the structures of (a)0- 
(ET) 2 X, (b)<9-(ET) 2 RbZn(SCN) 4 in the CO phase, (c)a-(ET) 2 I 3 
in the metallic phase, and (d) a-(ET)2l3 in the CO phase. The 
gray solid lines indicate the unit cell. For (b), (c), and (d), the 
transfer integrals estimated by the extended Hiickel method 11 ' 22 
are also shown. For a-(ET)2l3, the site indices A, B, C, and D 
in our notation correspond to A, B, C, and A f , respectively, in 
the conventional notation. 43 



unit cell in the c-direction and six transfer integrals ex- 
ist at low temperatures. In order to take account of the 
lattice effects, we consider three kinds of electron-lattice 



couplings, Si y 



s a , and that originate from the 



c- and a-axis molecular translations and molecular rota- 
tion, respectively. They have been introduced in the pre- 
vious papers, 40 ' 41 being suggested by the X-ray experi- 
ment that shows rather complicated displacements of ET 
molecules through the CO transition. 11 First, the c-axis 
translation alternates t c . This modulation gives t c \ and 
t C 2 in Fig. 1(b). Similarly, the a-axis translation induces 
the modulation of t p \ and t p %. Finally, the rotational de- 
grees of freedom are taken into account by the changes 
in t P 2 and t p ±. This type of modulation is expected to be 
important since the horizontal CO is formed on the t P 4 
chains with hole-rich molecules. In fact, the experimen- 
tal estimation 11 indicates that the dependence of transfer 
integrals on a relative angle (called elevation angle 11 ) of 
neighboring ET molecules is large and allows uniformly 
decreasing (increasing) \t P 2\ (\t P 4\) on the horizontally 
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where y c , y a , and y^ are modulations due to s c , s a , and 
S0, respectively. Here the signs in eq. (5) are chosen so 
that the transfer integrals realized experimentally at low 
temperatures correspond to y\ > for I = c, a and </>. 40 
In the calculations, we use t p = 0.1 (eV) and t c = —0.04 
in the high temperature phase for tij in eq. (1). Then, 
the modulations y^Vai and y^ are self-consistently deter- 
mined for each set of electron-lattice coupling constants. 

2.3 a-(ET) 2 h 

The transfer integrals for a-(ET)2l3 in the high- and 
low-temperature phases are shown in Figs. 1(c) and 1(d), 
respectively. The unit cell contains four molecules in 
both phases. According to the X-ray structural analy- 
sis, 22 sites A and D are equivalent to each other owing to 
the inversion symmetry in the metallic phase, while the 
symmetry breaks below the CO transition temperature. 
Sites A and B (C and D) become hole-rich (hole-poor) in 
the horizontal CO state. The large transfer integrals, t^i 
and £52? form a one-dimensional zigzag chain in the con- 
ducting layer. The extended Hiickel calculation indicates 
that they gradually increase with decreasing tempera- 
ture. 22 At the CO transition, tbi changes into t&\ and 
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Fig. 3. (Color online) (a) Modulation y c and (b)hole densities as 
a function of s c for U = 0.7, V c /U = 0.35, and V p /V c = 0.7. 



Fig. 4. (Color online) (a)Modulation and (b)hole densities as 
a function of s<t , for U = 0.7, V c /U = 0.35, and V p /V c = 0.7. 



£bi' in Fig. 1(d), while tb2 becomes £b2 and £b2'- The 
other transfer integrals with smaller values also change 
at the transition. The lattice modulation in this system 
is suggested to come mainly from the molecular rotation 
rather than the translational shifts of the molecules. 22 

In the present paper, we investigate the lattice effects 
on the system by taking account of the modulations in tbi 
and tb2- We do not consider the modulations in smaller 
transfer integrals than the above ones. These effects will 
be mentioned in the following section. As in the case 
of (9-(ET) 2 RbZn(SCN) 4 , we let s bl and s b2 denote the 
electron-lattice couplings that change tbi and ^2, respec- 
tively. The corresponding modulations are written as ybi 
and yb2. The transfer integrals in the low temperature 
phase are then written as 

tBv = hi + ybi , 

tBi = hi — ybi , 

(6) 

tB2> = U2 + Vb2 , 
tB2 = U2 — Vb2 , 

where the signs in eq. (6) are again chosen so as for pos- 
itive values to correspond to the modulations realized 
in a-(ET)2l3- The other transfer integrals are assumed 
to be unchanged. For the values of tij in eq. (1) with 
which the Hartree-Fock calculations are carried out, we 
use those in the high temperature phase that are shown 



in Fig. 1(c). 

3. Results and Discussions 

3.1 0-(ET) 2 RbZn(SCN) 4 

First, let us consider the effect of each electron-lattice 
coupling on the horizontal CO that is experimentally 
observed in <9-(ET) 2 RbZn(SCN) 4 . Figures 3, 4 and 5 
show the lattice displacements and order parameters as 
a function of the electron-lattice coupling when we as- 
sume the horizontal CO as the mean- field order parame- 
ter for U = 0.7, V c /U = 0.35, and V p /V c = 0.7. Note that 
the obtained CO state is not the ground state for small 
electron-lattice couplings, as we will discuss later. From 
Figs. 3(a) and 4(a), we can see that y c and increase 
with increasing s c and s^, respectively. The modulations 
grow linearly for small s c and s^, which indicates that the 
horizontal CO is stabilized by these lattice modulations. 
The hole density at each site in the unit cell shows that 
the disproportionation is enlarged by although the de- 
pendence is weak in the case of s c . On the other hand, 
the behavior of y a is different from the above two modu- 
lations, as shown in Fig. 5. y a becomes nonzero when s a 
exceeds some value. With the increase of y a , the order 
parameter tends to decrease, which shows the y a mod- 
ulation suppresses the charge disproportionation in the 
horizontal CO. The energy gain by these electron-lattice 
couplings is shown in Fig. 6. We find that the effect of 




Fig. 5. (Color online) (a) Modulation y a and (b)hole densities as 
a function of s a for U = 0.7, V c /U = 0.35, and V p /V c = 0.7. 
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Fig. 6. (Color online) Energy gain due to each electron-lattice 
coupling in the horizontal CO for U = 0.7, V c /U = 0.35, and 
V p /V c = 0.7. 



is the largest. Since the horizontal stripe CO is formed on 
the hole-rich sites that are connected by the t P 4 chain, the 
y<f> modulation lowers the energy through the exchange 
coupling between neighboring spins on the stripe. The 
energy gain due to spin fluctuations is also obtained by 
the perturbation theory from the strong-coupling limit, 



i.e., U 







Next, we discuss the stability of various CO patterns. 



Fig. 7. (Color online) (a) Ground-state energies for various CO 
states as a function of the ratio V p /V c for U = 0.7 and V c /U = 
0.35. The energy of 3-fold CO is chosen to be zero, hr, dg, and 
vt mean the horizontal, diagonal, and vertical COs, respectively, 
(b) Modulations in the transfer integrals for the horizontal and 
3-fold COs. 



The ground- state energies as a function of V p /V c are com- 
pared in Fig. 7(a). In this figure, the energy of the 3-fold 
CO without electron-lattice couplings is set at zero. We 
have shown only the lowest energy state for each CO pat- 
tern among different spin configurations. In the absence 
of electron-lattice couplings, the 3-fold CO with a ferri- 
magnetic spin configuration is the most favorable in the 
nearly isotropic region, V p /V c ~ 1, while the diagonal CO 
whose spin configuration is antiferromagnetic along the 
stripe and between stripes on the c-axis is stable when 
V p /V c is small, V p /V c < 0.7. 34 ' 40 For the horizontal CO, 
we plotted the energy of the state that is antiferromag- 
netic along the stripe and ferromagnetic between stripes 
on the c-axis. The state that is antiferromagnetic on the 
c-axis has a close energy and is nearly degenerate with 
the above state. Without electron-lattice coupling, the 
energy of the horizontal CO cannot be the lowest in any 
region. 

Figure 7(a) shows that the electron-lattice couplings 
stabilize the horizontal CO considerably. Here, we con- 
centrate on the effects of s c and separately since the 
results in the case where the three electron-lattice cou- 
plings are simultaneously present have been shown pre- 
viously. 40 The diagonal and vertical COs are not affected 
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by s c and s^, then the horizontal CO becomes more sta- 
ble than these states. The values of the electron-lattice 
couplings are chosen to be s c = 0.08 and = 0.1. The 
resultant lattice displacements are shown in Fig. 7(b), 
which are comparable to the experimentally observed 
values. The horizontal CO has hole-rich sites on the t P 4 
chains, which is consistent with the experimental find- 
ings. The 3-fold CO also has energy gain from the lattice 
modulations due to s c and s^. To obtain the mean-field 
solution for the 3-fold CO that coexists with lattice mod- 
ulations, we have assumed the unit cell of 2 x 6 sites in the 
a-c plane. In this state, a weak horizontal charge modula- 
tion is caused by the lattice distortion in the background 
of the 3-fold CO. The distortions in transfer integrals are 
weak compared to those in the horizontal CO. Note that 
the 3-fold CO is metallic even if the transfer integrals 
are modulated while the horizontal CO is insulating. In 
short, the horizontal CO with lattice distortion becomes 
stable for V p /V c < 0.7, while the 3-fold CO is favorable 
for 0.7 < V p /V c < 1. 

Here we briefly mention the effect of s a . Since the y a 
distortion reduces the horizontal CO, y c and tend to 
be suppressed when y a is large. 40 The energy of the hor- 
izontal CO is only slightly lowered by this y a distortion. 
The relative stability of this CO compared to the others 
in the range of V p /V c < 0.7 is unchanged as long as we 
choose s a to obtain the experimentally observed value of 

At finite temperatures, the stability of COs is investi- 
gated by calculating the free energy within the Hartree- 
Fock approximation. The used values of U, V c /U and 
the electron-lattice couplings are the same as those for 
T = 0. The temperature dependences of the free ener- 
gies for various CO patterns are shown in Figs. 8(a) and 
8(b) in the cases of V p /V c = 1.0 and V p /V c = 0.69, re- 
spectively. For V p /V c = 1.0, the 3-fold CO is the most 
stable state in a wide temperature range. On the other 
hand, the horizontal CO with lattice modulations be- 
comes more stable when V p /V c is small. In particular, a 
first-order metal-insulator transition between the 3-fold 
CO and the horizontal CO occurs for smaller values of 
V p /V c than unity, as shown in Fig. 8(b) where the tran- 
sition is located at T c ~ 0.025. 

The first-order metal-insulator transition at a finite 
temperature can be related to the experimental results 
of (9-(ET) 2 RbZn(SCN) 4 although the wave vector of the 
charge modulation in the 3-fold CO state is different from 
that of the experiments in the metallic phase. It has re- 
cently been proposed that the experimental observation 
can be reproduced by considering Coulomb interactions 
extended to a longer range than the nearest-neighbor 
sites. 36 As for the spin degrees of freedom, both the 3-fold 
and the horizontal COs in our Hartree-Fock calculation 
have spin orders which have not been observed in the 
experiments. The effect of quantum fluctuations is nec- 
essary in discussing the behavior of the spin degrees of 
freedom. 

The estimation of the intersite Coulomb interactions 
V p and V c indicates that these values are comparable, 44 
V p /V c ~ 1, where the 3-fold CO is the most stable in 
our calculation. A variational Monte Carlo study 35 for 
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Fig. 8. (Color online) Free energies for various CO patterns for 
(a) Vp/V c = 1.0 and (b)V p /V c = 0.69 in the case of U = 0.7 and 
V c /U = 0.35. pm means the paramagnetic metallic state with 
uniform charge density. A first-order metal-insulator transition 
from the 3-fold CO to the horizontal CO occurs at T c ~ 0.025. 



the extended Hubbard model without electron-lattice 
couplings also shows that the 3-fold CO is stable for 
Vp/Vc ~ 1. According to the recent exact-diagonalization 
study, 41 the horizontal CO with lattice distortions be- 
comes more stable even at V p /V c ~ 1 if we take ac- 
count of quantum fluctuations that are neglected in the 
Hartree-Fock approximation. It is also indicated that the 
Holstein-type electron-lattice coupling stabilizes the hor- 
izontal CO by an exact-diagonalization study. 31 

We find that the 3-fold state with coexisting weak hori- 
zontal charge modulation is stable at the nearly isotropic 
region V p /V c ~ 1. This can be related to the X-ray 
diffraction results on (9-(ET) 2 CsZn(SCN) 4 , 14 ' 15 which 
show two types of COs coexisting as short-range fluc- 
tuations. Although the present Hartree-Fock calculation 
gives an artificial long-range CO, quantum fluctuations 
are expected to destroy the long-range order and result 
in a state where different types of COs coexist as short- 
range fluctuations. 

3.2 a-(ET) 2 h 

As in the previous subsection, we first discuss the ef- 
fects of two electron-lattice couplings Sbi and Sb2 on 
the horizontal CO. In the CO state, sites A and B are 
hole-rich, while sites C and D are hole-poor. The spin 
configuration is antiferromagnetic along the stripe and 
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Fig. 9. (Color online) (a)Modulation y bl and (b)hole densities as Fi §- 10 - ( Color online ) (a)Modulation y b2 and (b)hole densities 
a function of s bl for U = 0.7, V c /U = 0.4 and V p /V c = 0.8. as a function of s b2 for C/ = 0.7, V c /U = 0.4 and y p /y c = 0.8. 



ferromagnetic between the stripes. We choose U = 0.7, 
V c /U = 0.4, and V p /V c = 0.8 here. Figures. 9 and 10 
show the lattice displacements and order parameters for 
the horizontal CO in the presence of s b i and 852, respec- 
tively. Although the charge distribution does not change 
largely, the displacements increase with increasing the 
electron-lattice coupling in both cases, which shows that 
these couplings stabilize the horizontal CO. This can be 
seen from the energy gain due to the distortions. As 
shown in Fig. 11, the energy of the horizontal CO is 
lowered in the presence of s b i and s b 2. 

The ground- state energies of various CO patterns per 
site are compared in Fig. 12(a). We set U = 0.7, V c /U = 
0.4 here. The electron-lattice couplings are chosen to be 
s bl = 0.08 and s b2 = 0.05. The 3-fold CO has a ferri- 
magnetic spin configuration in the 2x6 unit cell as in 
the case of the #-type salt. For the vertical CO, the en- 
ergies of two solutions that are nearly degenerate are 
shown, vtl means the vertical CO with hole-rich sites A 
and D whereas vt2 denotes the one with hole-rich sites 
B and C. For vtl, the spins on the stripes are antifer- 
romagnetic and those between the stripes are ferromag- 
netic. vt2 has a ferromagnetic spin configuration on the 
stripes and an antiferromagnetic one between the stripes. 
It can be seen from Fig. 12(a) that the 3-fold CO is the 
most stable near the isotropic repulsions, i.e., V p /V c ~ 1. 
When the anisotropy in V p /V c is large, i.e., V p /V c < 0.65, 
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Fig. 11. (Color online) Energy gain due to each electron-lattice 
coupling in the horizontal CO for U = 0.7, V c /U = 0.4, and 
V p /V c = 0.8. 



on the other hand, the horizontal CO is more stable 
than the other CO states even in the absence of electron- 
lattice couplings. This has been pointed out within the 
Hartree approximation. 29 This is in contrast to the 0- 
type salt where the horizontal CO does not become the 
ground state without electron-lattice couplings. There- 
fore, the band structure which originates from the vari- 



J. Phys. Soc. Jpn. 



> 



SB 



Full Paper 

(a) 



Author Name 



3fo!d hr - +- 

dg pm 

Vtl —A— vt2 — A — 

hr(s b i = 0.08. s b2 = 0,05) -o- 



-* + ' -jo-'' 
x ^- , J5 x x— 



' ,j0'- 



0.4 0.5 



0.6 



0.7 



0.8 



0.9 




0.005 



0.6 0.7 0.8 
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0.4. (b) Modulations in the transfer integrals for the horizontal 
CO. 



ety of transfer integrals in the high temperature phase 
is considered to have an important role to stabilize the 
CO in the a-type salt. In the presence of the electron- 
lattice couplings, the horizontal CO is more stabilized 
and it becomes the ground state for V p /V c < 0.7. No 
other state is affected by these electron-lattice couplings. 
In Fig. 12(b), we plotted the modulations in the trans- 
fer integrals, which are consistent with the experimental 
values, i.e., y^i ~ 0.02, and ~ 0.01 22 in the region 
where the horizontal CO becomes the ground state. In 
the parameter range we have chosen, the paramagnetic 
metallic solution has always a higher energy than the 
other CO patterns at T = 0. In this state, the charge 
distribution is not uniform due to the low symmetry of 
the crystal structure. 45,46 Site B is hole-rich and site C 
is hole-poor while the charge densities on sites A and D 
are the same in the metallic phase. This state is the so- 
called zero-gap semiconducting state where the energy 
band has two contact points at the Fermi energy in the 
Brillouin zone. 47,48 

In above calculations, we do not show the modula- 
tions in transfer integrals such as t a \ and £54. In fact, we 
find that the distortion in t a \ slightly lowers the energy 
of the horizontal CO although the distortion is smaller 
than those in t^i and t^- For the change in £54, the exper- 
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Fig. 13. (Color online) Temperature dependence of (a)free ener- 
gies for various CO patterns, (b) charge density at each site, and 
(c) modulations y^i and in the case of U = 0.7, V c /U = 0.4, 
and Vp/Vc = 0.6. The inset in (b) shows the behaviors of charge 
densities near the transition temperature for U = 0.7 (left) and 
U = 0.8 (right) with fixed values of V c /U = 0.4 and V p /V c = 0.6. 



imentally observed modulation does not become a self- 
consistent solution to the Hartree-Fock calculation. This 
is due to the fact that the distortion gains no energy since 
the increasing bond £b4' connects the hole-poor sites C 
and D. Therefore, we expect that the changes in tbi and 
tb2 have a major role in stabilizing the horizontal CO 
and the obtained results are qualitatively unchanged by 
other modulations. 

Finally, we consider the case of finite temperatures. 
As discussed above, the anisotropy in V p /V c determines 
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the stability of the CO states as in the case of the 6- 
type salt. For V p /V c ~ 1.0 the 3-fold state is stable, 
while the horizontal CO is stable when V p /V c is small. 
In the latter case, a phase transition occurs from the 
paramagnetic metallic state to the horizontal CO. The 
temperature dependence of the free energies of various 
CO patterns for V p /V c = 0.6 is shown in Fig. 13(a). 
We see that the transition takes place at T c ~ 0.076. 
The temperature dependences of the order parameters 
and lattice displacements are shown in Figs. 13(b) and 
13(c), respectively. In the paramagnetic metallic phase, 
the charge densities on sites A and D are the same and 
the modulations are absent, while these sites are distinct 
from each other in the horizontal CO phase accompanied 
by the modulations in the transfer integrals. If Holstein- 
type electron-lattice couplings were present, the charge 
disproportionation would make a difference among the 
site energies at sites A, B and C, even in the metallic 
phase. For T > T c , the contact points in the zero-gap 
state deviate from the Fermi energy, accompanied with 
thermally blurred small Fermi pockets. Below the transi- 
tion temperature, the spins on the stripes are weakly an- 
tiferromagnetic, and the spin order grows towards T = 0. 
The transition is of first order although the discontinu- 
ity in the order parameters is small. In the inset of Fig. 
13(b), we show the behaviors of hole-densities near T c 
for U = 0.7 and U = 0.8 by fixing the ratios V c /U = 0.4 
and Vp/Vc — 0.6. The first-order transition is more evi- 
dent in the case of large Coulomb interactions. Note that 
when we take much larger values of the electron-lattice 
couplings, T c becomes higher and the transition from the 
paramagnetic metal to the paramagnetic horizontal CO 
becomes continuous. We also find that, in the present 
case of U = 0.7 the CO transition becomes continuous 
without electron-lattice couplings, which indicates that 
the lattice effect can alter the order of the transition. 

A first-order transition from the charge disproportion- 
ated metallic state to the horizontal CO accompanied by 
the lattice distortion is consistent with the experimental 
findings in a-(ET) 2 l3. As for the spin degrees of free- 
dom, the antiferromagnetic spin order is obtained in the 
calculated CO state whereas the spin gap is experimen- 
tally observed below T c in the material. 17 Although this 
result is due to the Hartree-Fock approximation which 
does not take account of quantum fluctuations, the spin- 
singlet formation below T c can be expected by the fact 
that the spins on the horizontal stripes along £ B2 ' and £b3 
bonds form alternating Heisenberg chains as discussed by 
Seo. 29 Experimentally, a CO with long periodicity such 
as the 3-fold CO is not observed in this compound al- 
though the present calculation shows that it is stable for 
Vp/Vc ~ 1.0. As in the case of the #-type salt, the effects 
of fluctuations may further stabilize the horizontal CO 
even in the region of V p /V c ~ 1.0. 

4. Summary and Conclusions 

In the present paper, we study the role of the lattice de- 
grees of freedom to the formation of CO in the quasi- two- 
dimensional organic conductors #-(ET) 2 RbZn(SCN) 4 
and a-(ET) 2 l3. By taking account of Peierls-type 
electron-lattice couplings that cause the experimentally 



observed lattice modulation in each compound, we in- 
vestigate the relevant extended Hubbard model within 
the Hartree-Fock approximation. It is found that the 
electron-lattice couplings stabilize the horizontal CO for 
both compounds, which is consistent with the experimen- 
tal observations. For <9-(ET) 2 RbZn(SCN) 4 , the effect of 
the lattice distortion is crucial to realize the horizontal 
CO, since it does not become the ground state without 
electron-lattice couplings. The combined contributions of 
electron and lattice degrees of freedom to the CO transi- 
tion are in fact suggested by the recent dielectric permit- 
tivity measurement. 49 Our results indicate that the 
distortion, which corresponds to the molecular rotation, 
is the most important among the three couplings in order 
to obtain the horizontal CO. This result is also qualita- 
tively consistent with the exact-diagonalization study 41 
for eq. (1) on small clusters. At finite temperatures, a 
first-order metal-insulator transition is obtained, which 
is related to the CO transition in #-(ET) 2 RbZn(SCN) 4 . 
For a-(ET) 2 l3, on the other hand, the horizontal CO is 
already stable if we consider the full band structure and 
the anisotropy in the nearest neighbor Coulomb inter- 
actions, although the electron-lattice couplings further 
lower the energy of this CO. We have shown a finite- 
temperature first-order transition from the paramagnetic 
metallic state to the horizontal CO with modulations in 
the transfer integrals. This result agrees with the ex- 
perimental findings on a-(ET)2l3- Although the anti- 
ferromagnetic spin order is artificially obtained in the 
Hatree-Fock calculation, the spin gap behavior, which is 
observed in the material, is indeed expected since the ex- 
change couplings between the neighboring spins on the 
horizontal stripe are alternating as pointed out previ- 
ously. 29 

The contrastive role of the lattice degrees of freedom 
in each compound would appear in the different natures 
of the CO phase transitions. 6>-(ET) 2 RbZn(SCN) 4 in the 
metallic phase has a simple band structure with higher 
symmetry than that of a-(ET)2l3, and various COs com- 
pete with each other. When the 3-fold CO fluctuations 
are dominant at high temperatures as suggested by the 
present calculation, it is natural to expect that the tran- 
sition to the horizontal CO accompanies a large struc- 
tural distortion because it requires the considerable re- 
placement of the charge patterns with different unit cells. 
This results in the first-order transition with large discon- 
tinuity as observed experimentally. On the other hand, 
for a-(ET)2l3, the metallic state has already the charge 
disproportionation due to the complexity of transfer in- 
tegrals. Even in this metallic phase, site B is hole-rich 
while site C is hole-poor, which is common to the charge 
distribution in the horizontal CO. Because of this, the 
CO can be realized by merely breaking the equivalence 
of charge densities in sites A and D within the high- 
temperature unit cell. This can lead to the first-order 
transition with small discontinuity relative to that in the 
(9-(ET) 2 RbZn(SCN) 4 . 

Acknowledgment 

The authors would like to thank H. Seo and S. 
Miyashita for helpful discussions. This work was sup- 



10 J. Phys. Soc. Jpn. 



Full Paper 



Author Name 



ported by Grants-in-Aid for Scientific Research on Pri- 
ority Area "Molecular Conductors" (No. 15073224), for 
Scientific Research (C) (No. 19540381), for Creative Sci- 
entific Research (No. 15GS0216), and the Next Gener- 
ation Super Computing Project (Nanoscience Program) 
from the Ministry of Education, Culture, Sports, Science 
and Technology, Japan. 



T. Ishiguro, K. Yamaji and G. Saito: Organic Superconductors 
(Springer- Verlag, Berlin, 1998) 2nd ed. 

H. Seo, C. Hotta and H. Fukuyama: Chem. Rev. 104 (2004) 
5005. 

K. Miyagawa, A. Kawamoto and K. Kanoda: Phys. Rev. B 62 
(2000) 7679. 

R. Chiba, H. M. Yamamoto, T. Nakamura and T. Takahashi: 
J. Phys. Chem. Solids. 62 (2001) 389. 

Y. Takano, H. M. Yamamoto, K. Hiraki, T. Nakamura and T. 
Takahashi: J. Phys. Chem. Solids 62 (2001) 393. 
Y. Takano, K. Hiraki, H. M. Yamamoto, T. Nakamura and T. 
Takahashi: Synth. Met. 120 (2001) 1081. 

H. Mori, S. Tanaka and T. Mori: Phys. Rev. B 57 (1998) 12023. 
H. Tajima, S. Kyoden, H. Mori and S. Tanaka: Phys. Rev. B 
62 (2000) 9378. 

N. L. Wang, H. Mori, S. Tanaka, J. Dong, and B. P. Clayman: 
J. Phys.: Condens. Matter 13 (2001) 5463. 
K. Yamamoto, K. Yakushi, K. Miyagawa, K. Kanoda and A. 
Kawamoto: Phys. Rev. B 65 (2002) 085110. 
M. Watanabe, Y. Noda, Y. Nogami and H. Mori: J. Phys. Soc. 
Jpn. 73 (2004) 116. 

M. Watanabe, Y. Noda, Y. Nogami and H. Mori: J. Phys. Soc. 
Jpn. 74 (2005) 2011. 

R. Chiba, K. Hiraki, T. Takahashi, H. M. Yamamoto, and T. 
Nakamura: Phys. Rev. Lett. 93 (2004) 216405. 
M. Watanabe, Y. Nogami, K. Oshima, H. Mori and S. Tanaka: 
J. Phys. Soc. Jpn. 68 (1999) 2654. 

Y. Nogami, J. -P. Pouget, M. Watanabe, K. Oshima, H. Mori, 
S. Tanaka and T. Mori: Synth. Met. 103 (1999) 1911. 
K. Bender, K. Dietz, H. Endres, H. W. Helberg, I. Hennig, H. 
J. Keller, H. W. Schafer and D. Schweitzer: Mol. Cryst. Liq. 
Cryst. 107 (1984) 45. 

B. Rothamael, L. Forro, J. R. Cooper, J. S. Schilling, M. Weger, 
P. Bele, H. Brunner, D. Schweitzer and H. J. Keller: Phys. Rev. 
B 34 (1986) 704. 

H. Kino and H. Fukuyama: J. Phys. Soc. Jpn. 64 (1995) 1877. 
H. Kino and H. Fukuyama: J. Phys. Soc. Jpn. 65 (1996) 2158. 
R. Wojciechowski, K. Yamamoto, K. Yakushi, M. Inokuchi, 
and A. Kawamoto: Phys. Rev. B 67 (2003) 224105. 
21) S. Moroto, K.-I. Hiraki, Y. Takano, Y. Kubo, T. Takahashi, 
H. M. Yamamoto, and T. Nakamura: J. Phys. IV France 114 



22 
23 
24 
25 

26 

27; 

28 

29 
30 

31 

32 
33 

34 
35 

36 
37 

38 
39 

40 

41 

42 
43 

44 

45'' 

46 
47 
48 
49 



(2004) 399. 

T. Kakiuchi, Y. Wakabayashi, H. Sawa, T. Takahashi and T. 

Nakamura: J. Phys. Soc. Jpn. 76 (2007) 113702. 

T. J. Emge, P. C. W. Leung, M. A. Beno, H. H. Wang and J. 

M. Williams: Mol. Cryst. Liq. Cryst. 138 (1986) 393. 

H. Endres, J. Keller, R. Swietlik, D. Schweitzer, K. Angemund 

and C. Kruger: Z. Naturforsch. A 41 (1986) 1391. 

Y. Nogami, S. Kagoshima, T. Sugano and G. Saito: Synth. 

Met. 16 (1986) 367. 

Y. Nishio: unpublished. 

N. A. Fortune, K. Murata, M. Ishibashi, M. Tokumoto, N. 
Kinoshita and H. Anzai: Solid State Commun. 79 (1991) 265. 
S. Iwai, K. Yamamoto, A. Kashiwazaki, F. Hiramatsu, H. 
Nakaya, Y. Kawakami, K. Yakushi, H. Okamoto, H. Mori and 
Y. Nishio: Phys. Rev. Lett 98 (2007) 097402. 
H. Seo: J. Phys. Soc. Jpn. 69 (2000) 805. 

R. H. McKenzie, J. Merino, J. B. Marston, and O. P. Sushkov: 
Phys. Rev. B 64 (2001) 085109. 

R. T. Clay, S. Mazumdar and D. K. Campbell: J. Phys. Soc. 
Jpn. 71 (2002) 1816. 

T. Mori: J. Phys. Soc. Jpn. 72 (2003) 1469. 

J. Merino, H. Seo and M. Ogata: Phys. Rev. B 71 (2005) 
125111. 

M. Kaneko and M. Ogata: J. Phys. Soc. Jpn. 75 (2006) 014710. 
H. Watanabe and M. Ogata: J. Phys. Soc. Jpn. 75 (2006) 
063702. 

K. Kuroki: J. Phys. Soc. Jpn. 75 (2006) 114716. 

C. Hotta, N. Furukawa, A. Nakagawa and K. Kubo: J. Phys. 

Soc. Jpn. 75 (2006) 123704. 

C. Hotta and N. Furukawa: Phys. Rev. B 74 (2006) 193107. 
M. Udagawa and Y. Motome: Phys. Rev. Lett. 98 (2007) 
206405. 

Y. Tanaka and K. Yonemitsu: J. Phys. Soc. Jpn. 76 (2007) 
053708. 

S. Miyashita and K. Yonemitsu: Phys. Rev. B 75 (2007) 
245112. 

M. Watanabe: private communication. 

T. Mori, A. Kobayashi, Y. Sasaki, H. Kobayashi, G. Saito, and 

H. Inokuchi: Chem. Lett. 13 (1984) 957. 

T. Mori: Bull. Chem. Soc. Jpn. 73 (2000) 2243. 

A. Kobayashi, S. Katayama, K. Noguchi and Y. Suzumura: J. 

Phys. Soc. Jpn. 73 (2004) 3135. 

A. Kobayashi, S. Katayama, and Y. Suzumura: J. Phys. Soc. 
Jpn. 74 (2005) 2897. 

S. Katayama, A. Kobayashi, and Y. Suzumura: J. Phys. Soc. 
Jpn. 75 (2006) 054705. 

A. Kobayashi, S. Katayama, Y. Suzumura, and H. Fukuyama: 
J. Phys. Soc. Jpn. 76 (2007) 034711. 

F. Nad, P. Monceau, and H. M. Yamamoto: Phys. Rev. B 76 
(2007) 205101. 



